#
# $Id: scatter.dem,v 1.7 2003/10/28 05:35:54 sfeam Exp $
#
# Simple demo of scatter data conversion to grid data.
#

rgb(c) = (c==0) ? 65280 : (c==1) ? 255 : 16711680

set title "5S Knudsen-Hein Grammar Productions"
#unset hidden3d
set autoscale
#set style data points
set xlabel "(S->LS)"
set ylabel "(L->dFd')"
set zlabel "(F->dFd')"
set view 45,15  # Angle View 1
# set view 0,0  # Straight View 1

set term png
set output 'htmlImages/5S_khProductions1.png'

splot "dat/5S_khProductions.dat" using 1:3:5:(rgb($7)) with points pt 7 ps 1 lc rgb variable;


replot
set title "5S Canonical Ratio vs. Non-Canonical Ratio"
#unset hidden3d
set autoscale
#set style data points
set xlabel "Num Non-Canonical BPs/Seq. Length"
set ylabel "Num Canonical BPs/Seq. Length"

set term png
set output 'htmlImages/5S_2dBpRatio1.png'

plot "dat/5S_2dBpRatio.dat" using 1:2:(rgb($3)) with points pt 7 ps 1 lc rgb variable;


replot
set title "5S Watson-Crick vs. Wobble vs. Non-Canonical Ratio"
set autoscale
set xlabel "Num Non-Canonical BPs/Seq. Length"
set ylabel "Num W.C. BPs/Seq. Length"
set zlabel "Num Wobble BPs/Seq. Length"

set term png
set output 'htmlImages/5S_3DBpRatio1.png'

splot "dat/5S_3DBpRatio.dat" using 1:2:3:(rgb($4)) with points pt 7 ps 1 lc rgb variable;

replot
set title "5S Knudsen-Hein Grammar Productions"
#unset hidden3d
set autoscale
#set style data points
set xlabel "(S->LS)"
set ylabel "(L->dFd')"
set zlabel "(F->dFd')"
# set view 45,15  # Angle View 1
set view 0,0  # Straight View 1
fit m*x +c 'dat/5S_khProductions.dat' using 1:3:5 via m,c

set term png
set output 'htmlImages/5S_khProductions2.png'

splot "dat/5S_khProductions.dat" using 1:3:5:(rgb($7)) with points pt 7 ps 1 lc rgb variable;


replot
set title "5S Watson-Crick vs. Wobble vs. Non-Canonical Ratio"
set autoscale
set xlabel "Num Non-Canonical BPs/Seq. Length"
set ylabel "Num W.C. BPs/Seq. Length"
set zlabel "Num Wobble BPs/Seq. Length"

set term png
set output 'htmlImages/5S_3DBpRatio2.png'

splot "dat/5S_3DBpRatio.dat" using 1:2:3:(rgb($4)) with points pt 7 ps 1 lc rgb variable;


replot
set title "5S GC vs. AU vs. GU Basepair Ratio"
#unset hidden3d
set autoscale
#set style data points
set xlabel "p(dFd->gFc) / p(dFd->gFu)"
set ylabel "p(dFd->aFu) / p(dFd->gFu)"

set term png
set output 'htmlImages/5S_2DPfoldBpRatio.png'

plot "dat/5S_2DPfoldBpRatio.dat" using 1:2:(rgb($3)) with points pt 7 ps 1 lc rgb variable;


replot
set title "5S Nucleotide Distribution for U vs. C"
#unset hidden3d
set autoscale
#set style data points
set xlabel "p(s->U)"
set ylabel "p(s->C)"

set term png
set output 'htmlImages/5S_2DPfoldNucDist.png'

plot "dat/5S_2DPfoldNucDist.dat" using 1:2:(rgb($3)) with points pt 7 ps 1 lc rgb variable;


replot
set title "5S F-Measure Distribution"
#unset hidden3d
set autoscale
#set style data points
set xlabel "F-Measure"
set ylabel "Frequency"

set term png
set output 'htmlImages/5S_FMeasureDist.png'

plot "dat/5S_FMeasureDist.dat" using 1:2 with boxes fs solid 0.7 lc rgb "#050505"

replot
set title "5S F-Measure Distribution Easy"
#unset hidden3d
set autoscale
#set style data points
set xlabel "F-Measure"
set ylabel "Frequency"

set term png
set output 'htmlImages/5S_FMeasureDist_easy.png'

plot "dat/5S_FMeasureDist_easy.dat" using 1:2 with boxes lc rgb "green" fs solid 0.7

replot
set title "5S F-Measure Distribution Medium"
#unset hidden3d
set autoscale
#set style data points
set xlabel "F-Measure"
set ylabel "Frequency"

set term png
set output 'htmlImages/5S_FMeasureDist_medium.png'

plot "dat/5S_FMeasureDist_medium.dat" using 1:2 with boxes fs solid 0.7 lc rgb "blue"

replot
set title "5S F-Measure Distribution Hard"
#unset hidden3d
set autoscale
#set style data points
set xlabel "F-Measure"
set ylabel "Frequency"

set term png
set output 'htmlImages/5S_FMeasureDist_hard.png'

plot "dat/5S_FMeasureDist_hard.dat" using 1:2 with boxes fs solid 0.7 lc rgb "red"